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We show the stepped-pressure equilibria that are obtained from a generalization of Taylor relaxation known as 
multi-region, relaxed MHD (MRXMHD) are also generalizations of ideal MHD. We show this by proving that 
as the number of plasma regions becomes infinite, MRXMHD reduces to ideal MHD. Numerical convergence 
studies demonstrating this limit are presented. 



I. INTRODUCTION 

The construction of magnctohydrodynamic (MHD) 
equilibria in three-dimensional (3D) configurations is of 
fundamental importance for understanding toroidal mag- 
netically confined plasmas. The theory and numerical 
construction of 3D equilibria is complicated by the fact 
that toroidal magnetic fields without a continuous sym- 
metry are generally a fractal mix of islands, chaotic field 
lines, and magnetic flux surfaces. Hole, Hudson, and 
Dewar^ have proposed a mathematically rigorous model 
for 3D MHD equilibria that embraces this structure by 
abandoning the assumption of continuously nested flux 
surfaces usually made when applying ideal MHD. In- 
stead a finite number of flux surfaces are assumed to ex- 
ist in a partially-relaxed plasma system. This model, 
termed a multiple relaxed region MHD (MRXMHD) 
model, is based on a generalization of the Taylor re- 
laxation model 2 3 in which the total energy (field plus 
plasma) is minimized subject to a finite number of mag- 
netic flux, helicity and thermodynamic constraints. The 
model leads to a stepped pressure profile, with the pres- 
sure jumps across the barrier interfaces counterbalanced 
by corresponding jumps in the magnitude of the magnetic 
field. 

Although it might be expected that this MRXMHD 
model would reduce to ideal MHD in the limit of contin- 
uously nested flux surfaces, the discontinuous stepped- 
pressure profiles exhibited by this model make this unin- 
tuitive. In this paper we prove that the MRXMHD model 
does reduce to ideal MHD in the limit of continuously 
nested flux surfaces and provide supporting numerical 
evidence using the Stepped Pressure Equilibrium Code 
(SPECKS This demonstrates that the model proposed 
by Hole, Hudson, and Dewar - reduces to usual results 
such as ideal MHD in the integrable limit where contin- 
uously nested flux surfaces exist. 

In the next section we give a summary of the 
MRXMHD model and its solution for a finite number 
of plasma regions. In Section |III| we prove the main re- 
sult of the paper, that MRXMHD reduces to ideal MHD 
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FIG. 1. Schematic of magnetic geometry showing ideal MHD 
barriers Xi, and the relaxed plasma regions IZi. 



in the limit of continuously nested flux surfaces. This 
is followed by supporting numerical evidence examin- 
ing the convergence of SPEC to axisymmetric continu- 
ous pressure-profile solutions in Section [TV] The paper is 
concluded in Section fVl 



II. THE MULTIPLE-REGION RELAXED MHD MODEL 

As introduced previousljff^Hll the MRXMHD model 
consists of N nested plasma regions Hi separated by ideal 
MHD barriers li (see Fig. [T]) . Each plasma region is as- 
sumed to have undergone Taylor relaxation 2 -' to a mini- 
mum energy state subject to conserved fluxes and mag- 
netic helicity. The energy functional for the MRXMHD 
model can be written as 
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In each plasma region TZi the term Ui is the potential en- 
ergy, Si the plasma entropy, and Hi the gauge-invariant 
magnetic helicitySl. The plasma regions TZi are enclosed 
by flux surfaces, and are constrained to have helicity H® , 
plasma entropy Sf, poloidal flux Aip p ^ and toroidal flux 
Aipt,i- The Cp i and are circuits about the inner 
(<) and outer (>) boundaries of TZi in the poloidal and 
toroidal directions, respectively. 

Local minimum energy solutions to Q-Q are ob- 
tained by requiring the variation of W to be zero. With a 
fixed outer boundary In, these solutions have the forrrP^ 
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where ([5| applies in each plasma region TZi, ^ applies 
on each ideal interface li, n is a unit vector normal to the 
plasma interface li (see Figure [I]), and [[a;]] = a;,+i — Xi 
denotes the change in quantity x across the interface li. 



surfaces, and dH° and dS° are the corresponding con- 
straints. This energy functional is completed by expres- 
sions for the infinitesimal helicity dH and plasma entropy 
dS. 

The infinitesimal helicity dH follows from Q, 



dH = dr 3 A ■ B - di/) p f dl-A- dipt f dl-A 

Jc p (s) Jc t (s) 



(8) 



where C*(s) and C p (s) are toroidal and poloidal circuits 
along flux surface s. This may be further simplified by 
defining the enclosed flux functions 
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where ipt{s) and ip p (s) are the toroidal and poloidal fluxes 
enclosed by the flux surface s. 

Using (JsJ — ( 10 ) and the infinitesimal for dS with ^ 
gives the infinite-interface MRXMHD energy functional 
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III. THE CONTINUOUSLY NESTED FLUX-SURFACE 
LIMIT 

In this section we show that the MRXMHD model 
reduces to ideal MHD as the number of plasma re- 
gions increases. We begin by obtaining the limit of 
the MRXMHD energy functional (lJl-Q for continuously 
nested flux surfaces. 

We take the continuously nested flux surface limit of 
the MRXMHD energy functional @-@ by taking the 
limit as the number of plasma volumes N — > oo and 
the volume and enclosed fluxes of each plasma region 
approach zero. In this limit the MRXMHD energy func- 
tional becomes 
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where s is an arbitrary flux-surface label, dH and dS 
are the infinitesimal amounts of helicity and plasma en- 
tropy respectively between infinitesimally separated flux 



where H°(s) and S°(s) are the helicity and plasma en- 
tropy constraints. 

The variation of this energy functional is subject to the 
constraints 10 ) on the poloidal and toroidal fluxes 



enclosed by each magnetic surface. As discussed by Spies, 
Lortz, and Kaiser^, these constraints lead to the follow- 
ing relationship between the variation of the vector po- 
tential 5A and the variation of the interface positions 
<5x, 



n x SA = — (n • <5x) B, 



(12) 



where n is a unit vector normal to the flux surface. 

In the next section we first reproduce the result of 
Taylor^ demonstrating that in the absence of pressure 
the time- independent solutions of (11 1 are nonlinear Bel- 
trami fields. This result is then generalized to non-zero 
pressure in Section [HI B| 



A. Zero pressure limit 



The zero-pressure limit of (11) may be taken by set- 
ting p — ¥ 0, v(s) — > 0. In this limit, we need to consider 



3 



the variation of this functional with respect to the vec- 
tor potential, the positions of the flux surfaces, and the 
Lagrange multiplier /i(s). 

The variation 8fi(s) is independent of <5A and <5x, and 
may therefore be considered separately. Requiring the 
variation of W with respect to n(s) be zero enforces the 
hclicity constraint on each flux surface, 



SW\, = -\ J Sn(s) (dH-dH ) =0, 

or equivalently, H(s) — H°(s). 

The remaining variation of W with p = is 
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where s(x) is the flux surface label as a function of po- 
sition. The variation of the terms on the second line of 



(111 with fixed /i(s) is zero as i/jt(s) and ipp( s ) are given 
functions of the flux surface label s. 

The variation of s(x) can be obtained by defining 
s(x) = s(x) + 5s (x) to be the flux surface label after in- 
terface perturbation, and using the requirement that the 
perturbation doesn't change the label of a flux surface: 
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Using (17), the energy functional (11) may be written 
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where the perturbation of the magnetic field <5B has been 
written in terms of the perturbation of the vector poten- 
tial using SB = V x 5 A. 

This expression may be simplified using the relation 
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where Q is an arbitrary single-valued vector field, and 
Eq. ( 12 ) and the assumption that the outermost interface 



remain fixed (i.e. n • Sx = on the boundary) have been 
used. 



The relation ( 19 1 may now be used with Q-)B and 



Q — s- fiA to simplify (18), 
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The last line of (21 ) is zero: 



J dr 3 A ■ B(Sx ■ Vfi) = - J dr 3 A ■ (V/x X SA) , (22) 
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where (12) has been used, noting that V/i(s) || n. 
The variation SW has now been shown to be 

SW = J dr 3 SA -[VxB — fi(s)B] . 



(24) 



It is tempting to conclude from (24) that VxB 



/i(s)B, however this is not true in general. The flux con- 
servation condition ( 12 1 requires that 8 A ■ B = 0, hence 



SA is not a completely free variation. Requiring that 
the energy variation SW in (24) be zero for all possible 



variations only shows that the coefficients of independent 
variations are zero. 

The potential variation SA can be written in terms of 



independent variations using ( 12 1 



SA — Sx x B + nSA±, 



(25) 



where SA± is the remaining free variation of A, which is 
perpendicular to the flux surfaces. SA± is independent 
of Sx. 



Using (25), the energy functional variation (24) may 



be written as 



5W=J dr 3 [Sx ■ (-J x B) + 8A±n ■ J] 



(26) 



As Sx and SA± are independent, the time-independent 
solutions satisfy 



J x B = 0, 
n - J = 0. 



(27) 
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These two conditions imply that the current is parallel 
to the magnetic field 



V x B = A(x)B, 



(29) 



for some A(x). As the fields and currents are time- 
independent V ■ J = implies that B • VA = 0, hence A 
is constant on a field line. 



4 



Time-independent solutions of the infinite interface 
limit of the MRXMHD model without pressure are there- 
fore nonlinear Beltrami fields 



V x B = A(a)B, 



(30) 



where a labels the field line. This is the result of Taylor * 



One might have expected fi(s) to replace A(a) in (30) 
because for a finite number of interfaces the plasma in 
each volume satisfies V x B = /^B [see (Jjjb]. However, 
there are also surface currents on the interfaces between 
the plasma volumes. In the limit of an infinite number of 
continuously nested surfaces, the plasma volume current 
will have contributions both from the volume and surface 
currents of the finite- N case. Only if the surface currents 
in the finite- N case are zero should we expect X(a) to be 
equal to jtx(s). For example, the surface currents will be 
zero if the /i, are all equal, and in this case the solution 
is 



V x B = fiB, 



(31) 



with [i a constant. In this case \(a) is equal to /x(s). 

In the next section we consider the effect of pressure 
on the time-independent solutions of the infinite interface 
MRXMHD model. 



B. Non-zero pressure 

For non-zero pressure, the additional terms to the vari- 



ation (26 



that must be considered arc the variations of 



p terms in (111. 



the v anc 

The variation with respect to v{s) enforces the plasma 
entropy constraint 



SW\ U 
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or S(s) = S°(s). 

The variation with respect to pressure is 

SW\ = f *pUL.-„( s )V/7-i 
p J [7-1 7 
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As the variation Sp is independent of Sx and <5A, time- 
independent solutions satisfy 
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which implies that p is constant on a flux surface. 

The remaining additional term to the energy variation 
in the zero-pressure case is the variation of u(s) as the 
interface positions are varied. This term is 



(35) 



The gradient of v can be written in terms of the pressure 
p using pi], 
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The variation of the MRXMHD energy functional in- 
cluding pressure is ( 26 1 with the additional term ( 35 1 



SW=J dr 3 [Sx ■ (Vp - J x B) + 6A±n ■ J] 
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where (36 1 has been used. As the variations <5x and SA_[ 



are independent, the time-independent solutions of the 
infinite interface MRXMHD functional satisfy 



J x B = Vp, 
n J = 0, 



(38) 
(39) 



which are the equations for ideal MHD. In the limit of 
continuously nested flux surfaces, MRXMHD is equiv- 
alent to ideal MHD. In particular, in the axisymmetric 
N -> oo limit MRXMHD reduces to the Grad-Shafranov 
equation. 

In the next section we use SPEC 4 to illustrate the con- 
vergence of the MRXMHD model to axisymmetric ideal 
MHD with continuous pressure profiles. 



IV. NUMERICAL ILLUSTRATION 

The Stepped Pressure Equilibrium Code^ solves the 
MRXMHD model Q-Q for an arbitrary (finite) num- 
ber of plasma regions. We use this code to illustrate the 
results of the previous section by showing the numeri- 
cal convergence of SPEC to an axisymmetric continuous 
pressure-profile ideal MHD solution as computed by the 
VMEC codcP. 

The equilibrium is defined by a given, fixed outer 
boundary, the pressure and rotational-transform profiles 
as a function of the normalized toroidal flux, s — VVV'enci, 
where ip cnc \ is the total enclosed toroidal flux. For this 
comparison we choose ip cnc i = 2-7T in units where (1q = 1. 

For the numerical convergence study we choose the 
fixed outer boundary to be an axisymmetric torus with 
circular cross-section: 



fl = l.O + O.3cos(0), Z = O.3sin(0). 



(40) 



We define the equilibrium by choosing the pressure and 
rotational transform flux functions. The continuous pres- 
sure profile is selected to be 



p(s) =po(l-2 S + s 2 ), 



(41) 



where po is to be adjusted; e.g. po = for zero-beta. The 
continuous transform profile is selected to be 



= L a + (i e - L a )s, 



(42) 



where L a = (8 + 79)/ (9 + 7IO) and t e = (1 + jl)/(9 + 
7 10), and 7 = (1 + v / 5)/2 is the golden mean. This 
transform profile is selected to ensure that the rotational 
transform on the ideal interfaces in the MRXMHD model 
are noble irrationals. This ensures stability of the ideal 
interfaces^. 
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FIG. 2. Pressure profile used for demonstrating the N — > 
oo limit of MRXMHD. The continuous curve is the pressure 
profile in < |41[ ) which is used with VMEC, and the stepped 
profile is the N = 16 approximation used with SPEC. 

As input to SPEC, these profiles are discretized as fol- 
lows. For convergence studies as the number of plasma 
regions N — > oo, is it convenient to have the SPEC inter- 
faces equally spaced in s/s. So, for i = 1, . . . , N, we de- 
fine Si = \Ji/N and the interface transforms as Li = t(si). 
The pressure in each volume is constructed so that 

Pi f ds= f p(s) ds. (43) 

J Si— I J Si- i 

A discrete pressure profile, with N — 16 is shown in 
Fig. @ 

A comparison of the SPEC interfaces, for an W = 16, 
zero-/? case (i.e. pq = 0), is shown in Fig. [3j and for a 
high-/3 case in Fig. |4j For the high-pressure case, po was 
increased to give a Shafranov shift about one third the 
minor radius. 

To quantify the difference between the SPEC and 
VMEC solutions, we define a measure of the difference 
in geometry of a given magnetic surface as 

A = / <20|x V mec(0)-xspec(0)|, (44) 
Jo 

where x(#) is the intersection of the surface with the 
= plane, and 9 is the polar angle. 

Figure [5] shows A computed between the representative 
s = 1/4 SPEC interface and the corresponding VMEC 
interface as the number of interfaces N is increased up to 
the maximum afforded by computational limitations and 
expedience of N = 128. In particular, the convergence of 
the error is second order, A ~ A^ 2 , as shown in Fig. [5] 

V. CONCLUSION 

We have demonstrated that the Multiple-Region Re- 
laxed MHD model reduces to ideal MHD in the limit 
of an infinite number of plasma regions. In this limit, 
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FIG. 3. Comparison of magnetic surfaces: zero-pressure 
(po = 0); the SPEC interfaces are shown in thin lines (upper 
and lower half) and the VMEC surfaces are shown in thick 
lines (upper half only). 
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FIG. 4. Comparison of magnetic surfaces: high-pressure 
(po = 16); the SPEC interfaces are shown in thin lines (upper 
and lower half) and the VMEC surfaces are shown in thick 
lines (upper half only). 



the magnetic geometry is characterized by continuously 
nested flux surfaces. The appeal of MRXMHD is that for 
a finite number of plasma regions, only a finite number of 
flux surfaces are assumed to exist. The rest of the plasma 
may be characterized by smoothly nested flux surfaces, 
islands, chaotic fields, or some combination of these. In 
particular, the work of Hudson et al.^* demonstrates the 
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FIG. 5. Convergence: The error (A) between the continuous 
pressure (VMEC) and stepped pressure (SPEC) solutions are 
shown as a function of the number of plasma regions N for 
the s — 1/4 SPEC interface. The dotted line shows the zero- 
beta case (po = 0), and the solid line shows the high-beta case 
(po = 16). The grey line has a slope —2, the expected rate of 
convergence. 



application of SPEC to a DIIID equilibrium with a fully 
3D boundary in which magnetic islands form. In future 
work we will apply MRXMHD and SPEC to the RFX 
Quasi-Single Helicity state 12 in which two magnetic axes 
have been shown to form. 
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